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We study the influence of population imbalance on the pairing, spin, and density instabilities of 
a two component ideal Fermi gas after a sudden quench of interactions near a Feshbach resonance. 
Over a large region of parameters the pairing instability is dominated by finite momentum pairing, 
suggesting the possibility of observing FFLO-like states in the unstable initial dynamics. Long- 
wavelength density instabilities are found on the BCS side of the resonance, and are interpreted as 
a precursor of the phase separation expected at equilibrium. On the BEC side of the resonance, the 
pairing instability is present for scattering lengths that are larger than a critical value that is only 
weakly dependent on population imbalance and always smaller than the scattering length at which 
the Stoner-like spin instability occurs. 
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I. INTRODUCTION 

Ultracold atom systems host a web of phenomena that 
have made it possible to study quantum many body prob- 
lems in new ways that are often remarkably revealing. 
For example, studies of fermionic atoms have enabled 
a comprehensive investigation of the crossover between 
the BCS and BEC limits of the superfluid phase pQ. 
Crossover physics studies have been further enriched by 
exploring the fate of the superfluid phase as a function 
of the population imbalance between the species [5] and 
there is now a fairly complete understanding of the equi- 
librium phase diagram of this system [3] . 

One goal of ultracold fermion research is to shed light 
on the physics of interacting electrons. Atomic systems 
can often be prepared with less uncontrolled disorder 
and more control over the system parameters that deter- 
mine interaction strengths, facilitating comparisons be- 
tween theory and experiment. However the fact that in- 
teractions between electrons are normally dominated by 
Coulomb repulsion, whereas interactions between atoms 
are attractive at the densities of cold-atom systems, can 
stand in the way of these quantum simulation efforts. 

A particular challenge is the simulation of itinerant 
electronic magnetism, which is a combined consequence 
of strong repulsive interactions and Fermi statistics. The 
effective interaction strength between low-kinetic-energy, 
low-density atoms is strongly repulsive only when their 
attraction supports a shallow bound state. This obser- 
vation motivated an experiment |4J which explored the 
possibility of driving ferromagnetic instabilities [5j by us- 
ing a Feshbach resonance to suddenly increase the effec- 
tive repulsive interaction strength. However, the analogy 
with electrons is incomplete, mainly because of the pos- 
sibility to form weakly bound Feshbach molecules out of 
pairs of atoms. Indeed, recent experiments [B] and theo- 
retical studies [5] indicate that molecule formation is 
substantial even during the initial dynamics. 

Motivated by these studies, in the present theoretical 
work we study the initial linearized dynamics of density, 
spin, and pairing fluctuations of a system of atoms con- 



taining two fermionic species with unequal populations 
that is prepared in an ideal Fermi gas state and placed 
near a broad Feshbach resonance. We follow the ap- 
proach of Pekker et al. [7] , generalizing their work for the 
case of balanced populations to the case of imbalanced 
ones. 

We find that on the BCS side of the resonance the dom- 
inant instabilities are pairing and density phase separa- 
tion instabilities. We also observe that in the scattering- 
lcngth/spin-polarization space FFLO finite- wavevector- 
pairing fluctuations are more prominent than FFLO 
states are in the equilibrium phase diagram. This sug- 
gests the possibility of observing FFLO-like features in 
the initial dynamics of an ultracold gas after a quench. 
Our results for the density instabilities on the BCS side 
are in qualitative agreement with previous studies of un- 
stable collective modes in polarized systems |9]. On the 
BEC side we find that the critical scattering length for 
the appearance of a long-wavelength spin-density insta- 
bility remains larger than the critical scattering length 
to trigger the pairing instability. This finding indicates 
that a spin-density instability brought about by effec- 
tive repulsions, which attempts to phase separate spin 
species, if observed, will necessarily be accompanied by 
substantial pair binding into Feshbach molecules even in 
the polarized gas. 

Our paper is organized as follows. In Sec. [H] we discuss 
the initial pairing instabilities in a spin-polarized Fermi 
gas placed near a Feshbach resonance. In Sec. |III| wc 
turn to an analysis of longitudinal spin and density in- 
stabilities, which are coupled in spin-polarized systems. 



In Sec. IV we present a qualitative analysis of the in- 
stabilities of the transverse spin densities. Finally, we 
conclude in Sec. [V] with a brief summary and discussion 
of our findings. 



II. PAIRING INSTABILITIES 

Near a broad Feshbach resonance the vacuum interac- 
tion between low-energy fermions that are distinguished 
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by a quantum number (referred to below as spin) is con- 
trolled only by their mutual scattering length. As a con- 
sequence, any potential whose range is much smaller than 
the average interparticle separation can be used as an 



effective potential to describe the physics near the reso- 
nance, provided it is tuned so as to reproduce the phys- 
ical scattering length. The effective Hamiltonian of the 
system can be written in the form 



fc,CT crl,o-2 J 



(1) 



where £, a (k) = k 2 /2m — e F(T , and v(r — r') is an appro- 
priate pseudopotential. By performing a ladder sum, the 
scattering amplitude in the Fermi sea (the Cooperon) can 



be expressed in terms of the scattering amplitude in vac- 
uum and the occupation numbers of Fermi sea states [7] : 



C-\E, P) = £ (i + l ^ m (E + e Ft + e Fi )-^ ) + J 
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where E = cj^ + ui^ is the total energy measured with 
respect to the sum of the chemical potentials of the scat- 
tering pair and P = k^ + k^ is their total momentum. 
Particles of the same spin do not interact directly. 

Within time dependent mean-field theory it can be 
shown that the poles of the susceptibility associated with 
the pairing amplitude A P = Y,k( c \ p+k-\ c \p-ki> coin " 
cide with those of the Cooperon [TU] . The presence of 
poles in the upper half of the complex energy plane sig- 
nals pairing modes that grow exponentially during the 
initial dynamics after a sudden quench of the interaction 
potential. The equation C^ 1 (E,P) = thus determines 
the structure of unstable pairing modes, with the imag- 
inary part of the pole energy for a given momentum, 
lm(Ep), revealing the growth rate of the unstable mode 
in question. 

The implications of this criterion for pairing instabil- 
ities in initial dynamics are summarized in Fig. [T] We 
have parameterized the interaction strength by the di- 
mensionless product kpa, where we define kp for a po- 
larized gas as the Fermi momentum of an unpolarized 
gas with the same total density, i.e. 2k p = kp^ + kp^. 
In the latter expression kp-^n are the Fermi momenta of 
the two species of fermions. We also define a population 
imbalance parameter 

J= n 1 -n i= 4 t -4 l 
r H +n i k Ft + k Fl > 

Fig. [I] illustrates how the fastest growing pairing modes 
depend on these two parameters. 

For an unpolarized gas the pairing instability is present 
at all negative scattering lengths (on the BCS side of 



the resonance), and grows most rapidly for pairs with 
zero total momentum. Interestingly, the instability sur- 
vives over a finite region on the BEC side of the Fesh- 
bach resonance where the scattering length is positive. 
The absence of an instability for small positive scatter- 
ing lengths signals the metastability of the gas with re- 
spect to its pair-fluctuation channels. For a polarized 
gas we find that the instability can be dominated by fi- 
nite momentum pairing. This does not imply that the 
equilibrium state is a FFLO condensate of finite momen- 
tum pairs, but does indicate that the initial dynamics of 
pairing modes is dominated by finite momentum pairing. 

The large area of the region in the kpa-8 phase di- 
agram dominated by pairing at finite momentum con- 
trasts with the one associated with the FFLO phase in 
the equilibrium phase diagram (The FFLO state is be- 
lieved to occupy a very small region of the equilibrium 
phase diagram for a gas in three dimensions [3]). It falls 
roughly in the region of the equilibrium phase diagram 
in which phase separation between superfluid and excess 
majority-spin fermions is expected [3J. Therefore, if a po- 
larized system is prepared as a non-interacting gas on the 
BCS side far from the resonance, and the interactions are 
rapidly switched to the region where the pairing instabil- 
ity exists, our results suggest the possibility of realizing 
states with substantial finite momentum-pairing during 
the initial dynamics. An important caveat to this picture 
is that the pairing instability is expected to compete with 
density instabilities associated with phase separation, as 
we will discuss in the next section. 

When approached from the BEC side, the pairing in- 
stability onset is only weakly dependent on spin imbal- 
ance. The instability boundary can be interpreted as the 
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FIG. 1: (color online) (a) Pairing instability diagram. The 
red region in this diagram corresponds to those popula- 
tions imbalances and scattering lengths at which the dom- 
inant pairing mode is at finite momentum, the blue region 
is that dominated by zero momentum pairing and in the 
white regions there are no unstable pairing modes. The 
dashed-dotted line is the approximate threshold for the Fes- 
hbach molecule formation limited by the energy of holes left 
in the Fermi sea discussed in the text, which nearly coin- 
cides with the onset of the pairing instability on the BEC 
side, (b) Growth rate of the most unstable pairing mode 
Im(EM) as a function of the inverse scattering length across 
the resonance. The curves correspond to populations im- 
balances 8 = {0,0.2,0.4,0.6,0.8,0.9,0.96,0.99}. Pairing in- 
stabilities occur over a narrower Ufo, region at larger polar- 
izations, (c) Typical growth rates of the unstable pairing 
modes as a function of the pairing momentum, for S — 0.5 
and I/Zcfo = {— 0.4, — 0.3, 0.2}. Im(Ep) is the imaginary 
part of the Cooperon pole energy for a given pairing momen- 
tum, these curves illustrate how the instability evolves from 
being dominated by pairing at finite momentum at negative 
scattering lengths to zero momentum at positive scattering 
lengths. We conjecture that quenches of the scattering length 
that approach the resonance from the BCS side might bring 
into realization states with substantial finite momentum pair- 
ing during the initial unstable dynamics. 



threshold beyond which the energy gained by molecule 
formation can be absorbed by adding two holes to the 
Fermi sea [7j. Approximating the binding energy as 
~ 1 /ma 2 , the critical scattering length is found from 
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This threshold is shown in Fig. [T] as a dashed-dotted line, 
and nearly coincides with the one of the onset of pairing 



instability on the BCS side. In particular, the maximum 
energy of two holes created in an unpolarized Fermi sea 
is kp/m, whereas for the fully polarized Fermi sea it is 
kp^/2m = 2~ 1 / 3 k 2 7 /m. Thus the ratio of the critical 
scattering length at which pair formation occurs in an 
unpolarized gas to the corresponding critical scattering 
length in a fully polarized (FP) gas is only slightly less 
than one, 



app 



,-1/6 



0.89. 



(5) 



The actual ratio obtained directly from the Cooperon is 
do/o-FP ~ 0.75, with 1/kpCLo s» 1.07 and 1/kpapp rj 
0.80. When approached from the BCS side the pairing 
instability at full polarization does not occur until the 
Feshbach resonance is crossed. At full polarization, we 
can perform a similar estimate for the critical scatter- 
ing length. The pair binding in this case occurs at fi- 
nite momentum, between particles at the Fermi surface 
of the majority spin species and minority spin particles 
with zero momentum, so that the momentum of the pair 
is P = kp^. The energy of the pair is the sum of the 
binding energy and the center of mass kinetic energy, 
E = — 1/ma 2 + P 2 /Am, while the holes left in both Fermi 
seas have vanishing energy since they lie near the Fermi 
surfaces. It follows that 
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where we have referred to this critical scattering length as 
a Fp L ° because the pairing occurs at finite momentum. 
We have l/kpapp LO « 0.6, which deviates considerably 
from critical value calculated directly from the Cooperon 
l/kpapp LO w 0.4. This discrepancy originates from its 
proximity to the unitary regime, where the molecular dis- 
persion is greatly affected by the presence of the Fermi 
sea. 



III. DENSITY AND LONGITUDINAL 
SPIN-DENSITY INSTABILITIES 

We now turn to the instability analysis of density, 



Pg = E 
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and spin-density along the polar- 
ization axis, s zq = J2k,a c la a z a Ck+qa- The analysis is 
again done by looking for the occurrence of poles with 
positive imaginary parts in the corresponding suscepti- 
bilities. The responses of these quantities are coupled in 
the polarized gas and share the same unstable collective 
modes. We defer the analysis of the unstable collective 
modes of the transverse spin-densities to Sec. |IV| 

In time-dependent mean-field-theory, the spin and den- 
sity dynamics after a sudden interaction-strength quench 
have a random phase approximation (RPA) form. In this 
recipe, the bare interaction potential enters the "denom- 
inators" of the RPA susceptibilities. Often the short- 
range interactions in a dilute gas are modeled with a 
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FIG. 2: a) Diagrammatic representation of the integral equa- 
tion for the spin susceptibility vertex function for spin up. b) 
Relation between the spin vertices and the longitudinal spin 
susceptibility. 



contact pseudo-potential whose strength is proportional 
to the vacuum scattering length. This prescription, for 
which the RPA quenching picture is justified, would in- 
evitably fail near the resonance. Thus a conflict arises 
in the description of the unstable dynamics of spin and 
density near the resonance. With this on mind we will 
examine density and spin instabilities by assuming the 
RPA form for the dynamics, but with the bare inter- 
action replaced by the Cooperon. We believe that this 
approach, which has been previously employed to study 
spin instabilities in the unpolarized gas ^] , is a rea- 
sonable compromise, but we do not know of a systematic 
justification for it. 

The basic question we explore in this section is whether 
a free gas prepared far from the resonance can encounter 
a long-wavelength spin or density instability as the scat- 
tering length increases without reaching the critical scat- 
tering length that determines the boundary for the pair- 
ing instability. If a spin-density instability that caused an 
increase in the local spin-polarization occurred, it could 
inhibit pairing and allow the physics of ferromagnetism 
to be studied in a cold atom context. As discussed below, 
our calculations indicate that when the resonance is ap- 
proached from the BEC side the spin-density instabilities 
always occur within regions that are already unstable to 
pair formation. Even though we do not compare directly 
the growth rates of both instabilities, this finding is dis- 
couraging for the realization of an instability in which 
effective repulsive interactions would attempt to enhance 
locally the spin polarization of the gas. 

In a polarized gas the density and longitudinal spin 
responses are mutually coupled, but decoupled from the 
transverse spin response. The density and longitudinal 
spin susceptibilities can be expressed in terms of the spin- 
resolved density-density response functions 



functions as follows: 



Xaa{k)±Xsa{k) 



i I G^GAk + q^^k), (8) 
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where the +(— ) sign corresponds to the density (spin) 
vertex, q and k are four component vectors that include 
frequency and momenta components, G a {q) is the single- 
particle Green's function, and a denotes the spin oppo- 
site to a. The vertex functions satisfy the following self- 
consistent integral equations (see fig.[2|: 

T%>'(q,k) = lTi I G s (q')G 5 (q'+k)C{q+q'+k)Tr(q\k), 
Jq' 

(9) 

where C is the Cooperon. The poles of the vertex func- 
tions determine the unstable collective modes. 

Integral equations of the above type are familiar [12] 
from Fermi liquid theory. In the limit of small wavelength 
\k\ << kp the dominant contribution to the integral over 
q' in Eq.([9| is associated with low-energy Green's func- 
tion poles. This property can be exploited |12j by assum- 
ing that r is slowly varying near these poles. This allows 
to substitute 



G^G^k + q) 



C<r(q, k)S(k - e Fa )8(k - k Fcr ), 

(10) 

where ^ (q, k) is chosen so that its angular average re- 
produces the long-wavelength, low-energy behavior of a 
a free-particle (Lindhard) response function: 
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This approximation corresponds physically to the idea 
that the initial dynamics of the ideal gas is dominated 
by those states that are closest to it in energy, and that 
they in turn are formed from the gas state by making 
particle-hole excitations that live near the Fermi surface. 
The integral equation for the interaction vertices can then 
be converted into a matrix equation involving the expan- 
sions of r,£ and C into spherical harmonics. By truncat- 
ing the resulting equations at the s-wave term we obtain 
an expression for the vertices which resembles the RPA 
result, 
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where C s (uj) = hj d(k ■ k')C(uJ,kFfk + kp^k'), is the 
s-wave average of the Cooperon. We expect the s-wave 
truncation to be reliable for the stability boundary esti- 
mate at long wavelengths since in the limit of small (w, q) 
we have 



In particular, Xzz = Xtt _ X-l-t+X-l-l- _ Xn and X PP = Xn + 
X^T+Xi^+^Ti- The collective modes can be conveniently 
studied by defining the density, T£, and spin, r*, vertex 
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where v a = mkp a /2ir 2 is the density of states at the 
Fermi surface for the spin a. Assuming that the static 
susceptibility (u)/vp a q — > 0) correctly predicts the ex- 
istence of long-wavelength instabilities, the only non- 
vanishing component in the spherical harmonic expan- 
sion of £(q, k) in this limit is the s-wave term. It can 
then be argued that the instability boundary is correctly 
estimated by the s-wave truncation of Eq. ( 12 1. 



Our results for the density and longitudinal spin- 
density instabilities in polarized gases are summarized in 
Fig. [3] We find that, when the resonance is approached 
from the BEC side, the pairing instability always occurs 
for smaller scattering lengths than the spin instability. 
The spin-density instability occurs at smaller scattering 
lengths than those predicted when the Cooperon is re- 
placed by the contact pseudopotential, C — » Ana/m, 
valid in the dilute limit, kpa — > 0. When combined with 
RPA, this replacement would imply a Stoner instability 
at 1/kpa = 2/-7T for the unpolarized gas, whereas with 
the Cooperon we obtain 1/kpa s» 0.94 [7]. The differ- 
ence can be understood as being due to the regularized 
scattering, which produces stronger repulsions when the 
resonance is approached from the BEC side. Comparable 
enhancements were found in Monte Carlo studies [I3lll4j . 
Nevertheless, our results indicate that the population im- 
balance increases the critical repulsion required to trigger 
spin instabilities, and this critical repulsion remains al- 
ways larger than the corresponding one to produce pair- 
ing. We conclude that the preparation of a polarized ini- 
tial state does not lead to circumstances under which the 
Stoner instability of longitudinal spin modes is present, 
but the pairing instability is not. 

Even though the instability boundaries for the spin- 
density instability can be calculated consistently within 
our approach, the growth rates of these instabilities are 
reliable only near the critical parameters at which the 
instability appears. Inside the unstable region, beyond 
a critical scattering length, we have found discontinuous 
behavior of the unstable modes as a function of momen- 
tum and scattering length. Since Eq. ( 12 ) is justified 



only at small frequencies and wave- vectors, it is difficult 
to judge the extent to which such behavior is reliable. 
These discontinuities are also present in the density re- 
sponse of the unpolarized gas on the BCS side and are 
thus inherited by the longitudinal spin modes at finite po- 
larizations because of their coupling to the density. We 
will therefore focus only on the instability boundaries in 
this work. 

We first discuss the instability boundaries approached 
from the BCS side. As the resonance is approached at 
large polarizations, a long-wavelength density instability 
is always encountered at a lower scattering length than 
the corresponding critical scattering length for igniting 
pairing instabilities. At small polarizations (5 < 0.4) 
this boundary becomes increasingly close to the one of 
the pairing instability. The spin-density instability on 
the BCS side is brought about by effective attractive in- 
teractions. By examining the residue of the susceptibility 
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FIG. 3: (color online) (a) Growth rate of the density- 
longitudinal spin instabilities for 8 = 0.4 and scattering 
lengths l/k F a = {0.6, 0.65, 0.9}. (b) Growth rate of the 
density/spin instabilities for S = 0.8 and scattering lengths 
1/kpa = {—0.3, —0.25, 0}. (c) Instability phase dia- 
gram. The blue shaded region corresponds to the density- 
longitudinal spin instability region and for comparison the 
dashed lines depict the boundary of the pairing instability re- 
gion, the dotted line separates the region dominated by finite 
momentum pairing from that dominated by zero momentum 
pairing (see Fig. [T|. When the resonance is approached from 
the BEC side the spin-density instability always appears at 
a scattering length where the gas is unstable to pairing as 
discussed in the text. 



pole, we have verified that this instability corresponds to 
an unstable mode in which the majority spin and minor- 
ity spin densities grow with the same sign locally (that is, 
either they both decrease or both increase at any given 
point), with the amplitude of minority spin density larger 
than that of the majority spins. In other words, this in- 
stability attempts to create an alternation of regions with 
small population imbalance and regions with depleted 
minority spins. The relative amplitudes of the spin den- 
sities of the unstable mode near the instability boundary 
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(14) 



We view this as the precursor of the creation of regions 
where the pairing instability is locally developed, alter- 
nated with regions where the minority spins are depleted, 
explaining why the pairing and density instabilities are so 
close together. We expect that initial fluctuation bubbles 
with a growing density of paired particles should eventu- 
ally merge and phase separate from the excess majority 
spin particles. This instability is therefore compatible 
with the equilibrium state, which has phase separated 
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FIG. 4: (color online) (a) critical repulsion strength for the on- 
set of transverse spin instabilities (solid red line). The dashed 
line is the repulsion strength needed to sustain a spin imbal- 
ance in equilibrium (i.e. S/i — 0), which is finite at S = f. 
The dashed dotted line is the critical strength for the onset of 
longitudinal spin instabilities for the contact pseudopotential 
model. All lines converge to the Stoner criterion as 8 — > 0. (b) 
Typical growth rates of the unstable modes (imaginary part 
of the poles) of the transverse spin instability, (c) Real part of 
the collective modes dispersion. The unstable modes evolve 
continuously from the conventional spin waves in equilibrium 
as explained in the text and illustrated in Fig. [5] 



superfluid and normal gas regions, although our instabil- 
ity boundaries only indicate the appearance of linearly 
unstable modes and do not necessarily coincide with the 
boundaries between equilibrium phases. 

We note that at smaller polarizations, S < 0.4, the 
boundaries of finite momentum pairing and density in- 
stabilities almost merge. In this regime we conjecture 
that sudden quenches might bring into realization states 
with substantial finite momentum pairing in the initial 
dynamics. 

Several of our observations on the unstable dynam- 
ics on the BCS side are in qualitative agreement with a 
previous study of the unstable modes in the spinodal re- 
gion of the polarized system [9]. There, however, it was 
assumed that the dynamics of pairing was essentially in- 
stantaneous, so that the order parameter adjusted to its 
local equilibrium value after the quench, whereas the den- 
sity and spin modes would respond on slower time scales. 
Although this is enforced at arbitrarily large wavelengths 
due to conservation laws, the dynamics of pairing and 



density-longitudinal spin modes near the resonance oc- 
curs on similar time scales ~ e^ 1 at shorter wavelengths 
as illustrated in Figs. [T] and [3j 



IV. TRANSVERSE SPIN INSTABILITIES 

As discussed previously, near the resonance the mo- 
mentum and frequency dependence of the effective two- 
particle interaction is important. As will become clear in 
this section, the instabilities of the transverse spin den- 
sity occur at finite wavevectors, with the typical wave- 
length of the most unstable modes decreasing as the im- 
balance of populations increases. This situation compli- 
cates the analysis of the integral equations for the trans- 
verse spin vertices, since the Fermi liquid simplifications 
discussed in the previous are justified only at long wave- 
lengths and small frequencies. Nevertheless, important 
insights into the nature of the transverse spin instabil- 
ity can be gained with the use of a frequency-momentum 
independent pseudopotential V(r) = g5{r), and we will 
present them in this section. We leave open for future 
studies a more systematic treatment of this instability. 

The response of the x, y spin densities can be obtained 
in terms of the susceptibility X^t(9' w )' which describes 
the response of s^(q) = X)fc c l^ c fc+9t to a perturbation 
driven by its conjugate field. Within RPA [TSl [H] one 
obtains 



with 



d 3 k n^(k) — n_i(k + q) 



(15) 



(16) 



where, £ a (k) — k 2 /2m + gn 9 . The splitting between 
the dispersion of majority and minority spins is A = 
g(n^ — n±), and, correspondingly, the difference in chem- 
ical potentials reads as S/u, = /i-f — /xj, = Sep — A. 

Collective modes are found by solving 1 + gx^(q, w) = 
0. This same equation determines the dispersion of spin 
waves in a ferromagnet in equilibrium [15] . In Fig. [4] we 
depict the critical repulsion strength g cr it at which an 
instability in the transverse modes appears. Since no un- 
stable modes exist in equilibrium, the critical repulsion 
strength for the appearance of unstable collective modes 
for a given polarization must be larger than the interac- 
tion strength that is needed to sustain such polarization 
in equilibrium, i.e. g cr ;t > 3cq(^)- 

We find that the unstable modes have a finite momen- 
tum q > kp-f — kp^, and that they evolve continuously 
from the conventional spin-waves that the system would 
have in equilibrium. The instability exists whenever the 
spin wave spectrum, i? sw (g), coincides with minus the 
energy of a particle-hole excitation of minority into ma- 
jority spin, i.e. when E sw (q) = + q) - £).(*)), for 
some k such that k < kp^, kp^ < \k + q\. This can be 
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FIG. 5: (color online) (a) Equilibrium (i.e. Sfi — 0) spin 
wave dispersion and particle-hole continuum for S = 0.9. The 
boundary of the particle-hole excitation energies of majority 
into minority corresponds to the dashed line, the solid lines 
bound the negative of the energies of minority into majority 
excitations, and the blue line is the dispersion of spin waves, 
(b) Zoom in of Fig. (a). The spin waves cease to exist when 
they become degenerate with particle-hole excitations of ma- 
jority into minority spin, which allows them to decay into 
particle-hole excitations, (c) Spin wave dispersion at the crit- 
ical repulsion g cr it for the onset of transverse instabilities for 
S — 0.9. (d) and (e) correspondingly show the imaginary 
(growth rate of instability) and real parts of the spin wave 
dispersion at a repulsion strength gi/o = 1-25 for S — 0.9. 
The instability appears when the spin waves have minus the 
energy of minority-into-majority particle-hole excitations as 
explained in the text, and it would look like a spin-density 
wave of transverse magnetization with a spatial modulation 
corresponding to the inverse wavevector of the fastest growing 
mode ~ qm\x during the initial dynamics. 



seen as a necessary condition for the spontaneous excita- 
tion of magnons and particle-hole pairs while conserving 
total energy, momentum and spin. It is thus analogous 
to the critical interaction for pair formation discussed 
in the BEC side of the resonance, but instead of the 
particle-particle binding to form a molecule we consider 
particle- hole binding to form a magnon of spin — 1 which 
leaves a particle-hole excitation in the polarized Fermi 
sea of spin +1. The absence of unstable modes at long- 



wavelengths is thus a consequence of the Pauli exclusion 
principle which prohibits minority-into-majority excita- 
tions of momenta q < k F ^ ~ k F ^ [T7]- Fig. [5] illustrates 
how the spin wave spectrum evolves from being damped 
at equilibrium to unstable beyond the critical interaction 
strength. 



V. SUMMARY 

Our results imply that one should expect considerable 
differences in the initial unstable dynamics of polarized 
and unpolarized gases of ultracold fermions, and that fea- 
tures of the initial unstable state can be significantly dif- 
ferent from those expected from the equilibrium phase 
diagram in the population imbalance-scattering length 
coordinates. 

When the resonance is approached from the BCS 
side, finite momentum pairing instabilities and density- 
longitudinal spin instabilities are encountered. Pairing 
at finite momentum is the fastest growing pairing insta- 
bility over a wide range of parameters. This should be 
contrasted with the equilibrium FFLO phase region, con- 
jectured to occupy a rather small portion of the phase di- 
agram [3]. The density- longitudinal spin instability can 
be viewed as a precursor of the phase separation between 
a balanced superfluid and excess majority fermions, ex- 
pected to settle in at longer time scales. We predict 
an evolution of unstable modes as the population im- 
balance increases, with a transition between finite mo- 
mentum pairing at small polarizations (S < 0.4), and 
phase-separation modes at larger polarizations. 

When the resonance is approached from the BEC side, 
we find that the zero-momentum pairing instability al- 
ways occurs before the longitudinal spin-density insta- 
bility. The latter instability evolves from what would 
be the Stoner instability in the unpolarized gas. This 
finding seems discouraging for simulating a Stoner-like 
transition in a polarized gas of fermions. Nevertheless, 
as illustrated with a simplified model in Sec. |IV[ insta- 
bilities of the transverse spin modes are also expected to 
appear before the longitudinal spin and density instabil- 
ities. Therefore, an analysis of these instabilities which 
systematically accounts for the modification of the in- 
teractions near the resonance is necessary to resolve the 
fate of the competition between the transverse-spin and 
pairing instabilities at finite polarizations. 

It is worth emphasizing that in a polarized gas with 
globally conserved populations of each spin species, some 
of the physics of itinerant electron ferromagnetism can be 
explored, independently of the presence of spin-density 
instabilities. For example, in our simple model of Sec. |IV[ 
the spectrum of spin waves acquires a positive mass 
(magnetization stiffness) for interaction strengths below 
any spin instability, and below the interaction strength 
that would spontaneously sustain the population imbal- 
ance in equilibrium. Experimental studies of spin-wave 
dynamics of fully polarized states close to a Feshbach res- 
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onancc are likely to shed light on many-electron physics 
questions that arise in the context of itinerant electron 
magnetism, and also on pairing physics questions that 
are unique to the cold atom problem. 
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